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£N) ■ The recent radial velocity Keck data for GJ876 in Laughlin et al. 

(2005) is shown to be in good agreement, apart from a 6 m/s scat- 
' ter, with a theoretical calculation (Nauenberg, 2002) based on orbital 

parameters from a fit to the earlier Keck data. The time variation of 
the periods of the inner and outer planets, which are locked in a near 
Q-f 2:1 resonance, are evaluated, and their mean values, Pj and P Q , are 

O ' shown to satisfy closely the relation P /Pi = 2 + P /P, where P is the 

period for the retrograde precesion of the commom mean periastron 
of each of these planets. 



Recently, Laughlin et al. (2005) have given a progress report on the 
extrasolar planetary system GJ976, which consists of two Jupiter size planets 
in a near 2:1 resonance, and provided new high precision radial velocities 
obtained with the Keck telescope during the period 2001 to 2004. This new 
report discusses a fit to a coplanar configuration of the two planets which 
includes a 6 meters/second stellar jitter. The possible occurrence of such 
a jitter was suggested in (Nauenberg 2002) because of the excesively large 
values of the reduced chi-square in the fit to the data, 1 , and it appears 
to have further confirmation by studies of chromospheric activity in M3-M4 
stars similar to GJ876 (Laughlin et al. 2005). There also may be, however, 
other sources, e.g. additional smaller planets of order 6/200 of a Jupiter 

1 In earlier fits to the data ( Marcy et al. 2001; Laughlin & Chamber 2001; Rivera & 
Lissauer 2001), the reduced chi square was reported incorrectly 
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mass which can also contribute to the estimated 4 to 6 m/s variations in the 
maximum radial velocity of about 200 m/s 2 . In Figs. 1-3 we show that the 
radial velocities in the new Keck data are in good agreement, apart from a 
mean 6 m/s scatter, with the values predicted from theoretical calculations 
(Nauenberg 2002) based on the orbital parameters obtained from a fit to the 
early Keck data of Marcy et al. (2001) including an additional 9 data points 
released previously. The scatter is shown in Fig. 4. This good agreement is 
surprising in view of the fact that in the new report (Laughlin et al. 2005) 
the velocities from the earlier data have been modified. This modification 
includes a substantial shift in the mean velocity of the central star of about 
70 meters/second which is included in comparing the new data with our 
predictions. 

It has been customary to characterize the orbital motion of the two plan- 
ets of GJ876 by the parameters of the osculating ellipse at a given epoch. 
In particular, the two parameteres Pj and P Q associated with the periods of 
the inner and outher planets are approximately 30 and 60 days, respectively, 
which was the original basis for the assertion that these two planets move in a 
near 2:1 resonance. But the actual periods, defined as the time for a sidereal 
revolution of the planets, are not constants, but exhibit oscillations 3 . The 
inner planet moves alternatively between two nearby orbits, Fig. 5, with two 
different oscillating periods which are shown in Fig. 6 . The oscillations of 
the period of the outer planet is somewhat more regular and is shown in Fig. 
7. These oscillations have a period of about 660 days. Of special interest is 
the mean value of the orbital periods, because the condition of a near 2:1 
resonance implies the relation 

P P 

where P is the mean period for the rotation of the periastron of these planets, 
which is moving retrograde relative to the rotation of the planets . Since this 
relation, which reveals the reason why P Q /Pi > 2, has not been presented 
previously, a derivation is given in the following Appendix. We have calcu- 
lated the mean values of these periods by counting the number of rotations 

2 After the completion of this paper, a new planet of 7.5 earth-mass orbiting GJ876 has 
been reported (Rivera et al. , to be published) 

3 The oscillations of the period, defined as the time elapsed between consecutive peri- 
astrons of the orbit, is shown (Nauenberg 2002) 
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for a given time, and increasing this time t until there are no further changes 
in the period to a chosen accuracy. Increasing t from 1000 to 8000 days, we 
obtain P Q = 61.035 days, Pj = 30.220 days and P = 3125 days. Hence, we 
obtain P D /Pj — 2 = .0197 and P /P = .0195, which is in very good agreement 
with our period relation, Eq. 1. Previously we had estimated that P = 3200 
days (Nauenberg 2002), while Lee and Peale (2002) obtained P = 3100 days, 
and more recently Laughlin et al. (2005) reported P = 3205 days (quoted 
as a rotation rate of 41° per year). Actually, we have evaluated P from the 
orbital parameters given in (Laughlin 2005), and obtained instead P = 3136 
in better agreement with our results. 

This analysis indicates that GJ876 offers a remarkable opportunity for 
a detail quantitative study of a near 2:1 resonance in a planetary system, 
which will greatly improve as more data becomes available in the future. 



Appendix, Period relation for a near r : s resonance 

In this appendix we consider the general condition that two planets are in a 
near r : s resonance, where r and s are two integers. Let u>i(t) and u> (t) be 
the time dependent angular frequencies for the inner and outer planets, and 
Qi(t) and Q {t) the corresponding angular frequencies for the rotation of the 
periastron of these two planets. Then the number of rotations n-i and n Q of 
the inner and outer planet from one periastron 4 to the the next one, during 
a given time t, is given by 

m = {uo t - Sli)t/2n, (2) 

and 

n D = (uj - Q )t/27i (3) 

where u>i, u , fli and Q Q are the mean values of these angular frequencies 
during t, 

Ui = \ f dt'cjiit'), (4) 

t Jo 

ooo = - f dt'u (t'), (5) 

t Jo 

4 The periastron is the position on the orbit nearest to the central star. It is assumed 
that this point is connected to the farthest point on the orbit by a line, referred to as the 
line of apsides (Lee 2000; Laughlin 2005). However, when the mean periastron is rotating, 
as in GJ876, these two points are found to be 2° to 3° from such an imaginary line 
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sii = l f dt'n t (t'), (6) 

t Jo 

and 

n a = j f dt'n (t'). (7) 

t Jo 

The condition for a near r:s resonance leads to an asymptotic condition for 
large t, which requires that 

n - = ~, (8) 

n Q s 

and 

Q i = n o = Q. (9) 

Hence, according to Eqs. 2 and 3, Q = ru Q — siVi)/(r — s) which for a 2:1 
resonance gives f2 = 2uo Q — uji (Lee 2002). For negative values of fi, which 
corresponds to retrograde motion for the mean position of the periastron of 
both planets, the mean periods Pi = 2n/ui, P Q = 2tt/uj and P = 2n/\fl\ 
satisfy the relation 

P fr r -s + ^T- < 10 > 

For the case of GJ876, r = 2 and s = 1, we obtain Eq. 1 

The angular oscillations of the periastron of the inner and outer planets 
relative to its mean precession, which have been discussed in (Lee 2002;Laugh- 
lin 2005), is given by 

9 t = C dt'iVtiit') - Vt) (11) 

Jo 

e = f dt'(n (t') - n) (12) 

Jo 

Our numerical results are shown in Figs. 8 and 9, giving maximum dis- 
placements \9i\ — 6° and \6 \ = 50° which are somewhat different from those 
reported in (Laughlin 2005). 
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Figure 1: New Keck radial velocities (Laughlin et. al 2005) and calculation 
based on a fit (Nauenberg 2002) to the early Keck data (Marcy et al. 2001) 
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Figure 2: New Keck radial velocities (Laughlin et. al 2005) and calculation 
based on a fit (Nauenberg 2002) to the early Keck data of Marcy et al. (2001) 
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Figure 3: New Keck radial velocities (Laughlin et. al 2005) and calculations 
based on a fit (Nauenberg 2002) to the early Keck data Marcy et al. (2001) 
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Figure 4: Scatter of new Keck radial velocities from a fit (Nauenberg 2002) 
to the early Keck data of Marcy et al.(2001) 
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Figure 5: Two consecutive orbits of the inner planet (red curve), and during 
a single orbit of the outer planet (green curve). The black lines indicate the 
positions of periastron for the inner planet at J, and the outer planet at O 
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Figure 6: Oscillations of the two periods of the inner planet in its consecutive 
orbits, indicated by triangles and squares. 
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Figure 7: Oscillations of the period of the outer planet 
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Figure 8: The periastron angle for the two consecutive orbits of the inner 
planet shown by triangles and squares 
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Figure 9: The periastron angle for the outer planet 
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